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ABSTRACT 

We examine the radial entropy distribution and its scaling using 3 1 nearby galaxy clusters from the Representative XMM-Newton 
Cluster Structure Survey (REXCESS), a sample in the temperature range 2-9 keV selected in X-ray luminosity only, with no bias 
toward any particular morphological type. The entropy profiles are robustly measured at least out to i^iooo in all systems and out to 
i?500 in thirteen systems. Compared to theoretical expectations from non-radiative cosmological simulations, the observed distributions 
show a radial and mass-dependent excess entropy, such that the excess is greater and extends to larger radii in lower mass systems. 
At Rsm, the mass dependence and entropy excess are both negligible within the large observational and theoretical uncertainties. 
Mirroring this behaviour, the scaling of gas entropy is shallower than self-similar in the inner regions, but steepens with radius, 
becoming consistent with self-similar at R^oq. There is a large dispersion in scaled entropy in the inner regions, apparently linked 
to the presence of cool cores and dynamical activity; at larger radii the dispersion decreases by approximately a factor of two to 
30 per cent, and the dichotomy between subsamples disappears. There are two peaks in the distribution of both inner slope and, 
after parameterising the profiles with a power law plus constant model, in central entropy Ko. However, we are unable to distinguish 
between a bimodal or a left-skewed distribution of Ko with the present data. The distribution of outer slopes is unimodal with a median 
value of 0.98, and there is a clear correlation of outer slope with temperature. Renormalising the dimensionless entropy profiles by 
the gas mass fraction profile /gas(< R), leads to a remarkable reduction in the scatter, implying that gas mass fraction variations with 
radius and mass are the cause of the observed entropy structural and scaling properties. The results are consistent with the picture 
of a cluster population in which entropy modification is centrally concentrated and extends to larger radii at lower mass, leading 
to both a radial and a mass-dependence in the gas mass fraction, but which is increasingly self-similar at large radius. The observed 
normalisation, however, would suggest entropy modification at least up to Rmo, and even beyond, in all but the most massive systems. 
We discuss a tentative scenario to explain the observed behaviour of the entropy and gas mass fraction in the REXCESS sample, in 
which a combination of extra heating and merger mixing maintains an elevated central entropy level in the majority of the population, 
and a smaller fraction of systems is able to develop a cool core. 

Key words. Cosmology: observations. Cosmology: dark matter, Galaxies: cluster: general, (Galaxies) Intergalactic medium. X-rays: 
galaxies: clusters 



1. Introduction 

The first order model of structure formation - that of hierarchi- 
cal, dark matter dominated gravitational collapse - is capable of 
reproducing only the gross statistical properties of the galaxy 
cluster population. In this scenario, the intracluster medium 
(ICM) is heated to X-ray emitting temperatures by shocks and 
compression as it falls into the potential well of the dark mat- 
ter, and the resulting X-ray cluster population is self-similar and 
scale-free. In real clusters, second order effects, linked primarily 
to feedback from galaxy formation and radiative cooling of the 
gas, serve to modify the X-ray p roperties of the ICM with respect 
to the se expectations (see, e.g.. IVoitll2005l : iBorgani & Kravtsovl 
l2009l for recent reviews). The effect of these nongravitational 
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processes is substantial. It can be seen most readily in the re- 
lation between the X-ray luminosity and temperature, which in 
the first order scenario scales simply as L 25 T^, but which is ob- 
served to scale as L'^T^ (e.g. Pratt et al.ll2009l and references 
therein), implying a progressive suppression of luminosity in 
low temperature systems. 

In recent years, spatially resolved observations have allowed 
us to examine in more detail the impact of nongravitational pro- 
cesses on the ICM, mainly through radialprofiles and mapping. 
In this context, the entropy K of the IQVy is of considerable in- 
terest because the observable X-ray characteristics of a cluster 

' Keeping with convention, we use the X-ray astronomer's 'entropy' 
throughout this paper Defined as A" = kT/nl^^, where n^ is the electron 
number density, this quantity is related to the true entropy by a loga- 
rithm and an additive constant. 
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Table 1. Basic cluster data. K(Rs) is the entropy measured at the radius corresponding to density contrast S. Cool core and morphological 
classification from lPratt et alj ( 120091) . 
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" Temperature in keV, measured in the [0.15 - 0.75] R^qo aperture 

* Mass in units of h^g 10'* Mq, estimated from the M500 - i^x relation given in Egn.fTl jArnaud et al.ll2009l) . 



are just manifestations of its distribution in the dark matter po- 
tential well. Entropy is generated during the hierarchical assem- 
bly process, yet is modified by any other process that can change 
the physical characteristics of the gas. It is thus a quantity that 
preserves a record both of the accretion history of a cluster and 
of the influence of non-gravitational processes on the properties 
of its ICM, and as such it is a useful tool for our understanding 
of the thermodynamic history of the cluster population. 

Early measurements of the entropy based primarily on 
ROSAT and ASCA data indicated that groups had flatte r en- 
tropy profiles than cluster scale objects ( David et al.lll996 ^. and 
measurements of the entropy at 0. 1 /?200 revealed an entropy- 
temperature (K - T) relation that was shaUower than expected 
jPonman et alJ 119991; iLlovd-Davies et alJ I2OOOI; iPonman et all 
120031) . These data also afforded the first indications for excess 
entropy above that expected from gravitational collapse even 
at large radius in g roup-scale objects (Finoguenov et al. 2002; 
IPonnian et al.ll2003l) . The advent of XMM-Newton and Chandra 
has allowed relatively high resolution spatially resolved mea- 
surement of the entropy t o be obtained across a wide range of 
cluster and group masses (Pratt & Arnaud 2003; Piffaretti et al.' 
2005; Pr att et al. 2006; Morandi & Ettori 2007; Na gai et al 
2007; Z hang etalJ 120081; ISanderson et al J 120091: IJohnson et al. 
2009i) . Recent results have suggested that the entropy is indeed 



highe r than expected from gravitational collapse at least out to 
R25Q0 ("Pratt et al."2006- Nagai et al. '2007; 'Sun et al."2009), and 
perhaps further (Pratt etal. 2006; Sun et al. 2009), even up to 
relatively high masses. In addition, indications for excess en- 
tropy h ave been found at l arge radius in intermediate redshift 
groups dJeltema et al.ll2006h . and the fir st measurements o f the 
entropy evolution have been undertaken dEttori et aLll2004l) . 

In the present paper we re-i nvestigate the entropy with 
R EXCESS dBohringer et alj l2007h . a representative sample of 
33 local (z < 0.2 ) clusters drawn from the REFLEX catalogue 
dBohringer et al J 12004,) . all of which have been observed with 
XMM-Newton. The properties of the REXCESS sample allow us 
to define a robust local reference for entropy structure and scal- 
ing. REXCESS was designed to be representative of any high- 
quality local X-ray survey, thus clusters have been selected in 
luminosity only, ensuring no morphological bias, in such a way 
as to sample the X-ray cluster luminosity function in an optimal 
manner. Moreover, distances were optimised so that R^oo falls 
well within the XMM-Newton field of view, increasing the preci- 
sion of measurements at large radii as compared to more nearby 
clusters, which often fill the field of view and for which back- 
ground modelling is consequently more complicated. 

In the following, we first examine the normalisation of the 
entropy with respect to predictions from non-radiative cosmo- 
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logical simulations - such 'adiabatic' simulations include only 
gravitational processes - finding a systematic entropy excess that 
is greater at small radii and in lower mass systems. The consider- 
able dispersion at small radii appears linked to whether a cluster 
possesses a cool core or is morphologically disturbed. The mass 
dependence disappears at Rsoo, implying that entropy scaling is 
self-similar, with a normalisation that is approximately consis- 
tent with predictions. Parameterising the profiles in terms of a 
power law plus constant model, there are two peaks in the distri- 
bution of central entropy but there is no strong evidence that it 
is bimodal. The distribution of outer slopes is unimodal and the 
slope depends on temperature. We then link the entropy scaling 
and structural properties to a systematic variation in gas content 
with total mass and with radius. Finally, we discuss mechanisms 
which could bring about the observed entropy characteristics, 
and propose a tentative scenario to explain the observed entropy 
distributions. 

We adopt a ACDM cosmology with Hq = 70 km s"' Mpc"' 
(i.e., hjo =1), Q.M - 0.3 and Qa = 0.7. All uncertainties are 
quoted at the 68 per cent confidence level. 

2. Sample and analysis 

2.1. Sample description and subsample definition 

A full description of the R EXCESS sample, including 
the XMM-New t on ob servation details, can be found in 
iBohringer et al.l (l2007h. and t he pre hminary data analysis 
is described in ICroston et al.l ( l2008l) . Two of the objects, 
RCXC J0956.4-1004 and RXC J2157.4-0747 (the Abell 901/902 
supercluster and a bimodal cluster, respectively), display com- 
plex morphologies which preclude their use for the present radial 
profile analysis. Basic cluster parameters are listed in Table[T] 

On occasion in the following, we will subdivide the sam- 
ple into cool core and non-cool core systems, or according to 
whether the clusters are morphologically relaxed or unrelaxed. 
These subsamples were established to cull approximately the 
most extreme thirty perc ent of the full sample in each cate- 
gory and are defined as in lPratt et al.l (2009). Thus clusters with 
central density EizT^iiffi > 4 x 10"^ cm""* are classifiecQ as 
cool core systems (10/31), and those with centre shift parameter 
(w) > 0.01 /?5oo (derived with the central regions excised) are 
classified as morphologically disturbed (12/31). Both the cen- 
tral densities w^.o and centroid shift parameter (w) are given in 
iHaarsma et al] (l2009l) . 

2.2. Data analysis 

Event lists were processed using version 7.0 of the 
XMM-Newton SAS. All data products were extracted from 
event lists that were generated, cleaned, PATTERN-selected, 
vignetting-correc ted, and point source-removed as described in 
iPratt et all (|2007|) . 

2.2.1 . Gas density profiles 

The procedure used to calculate the gas density profiles, plus 
extensive ana l ysis o f their properties, is described in full in 
ICroston et al.l (l2008h . In brief, surface brightness profiles, cen- 
tred on the peak of the X-ray emission, were extracted from 3'.'3 
bins in the [0.3-2] keV band and deprojected and PSF corrected 

^ E{z) is the ratio of the Hubble constant at redshift z to its present 
value, E^{z) = DmCI + z)' + Ha- 



using the non-parametric method introduced in ICroston et al.l 
( l2006h . These were converted to gas density by calculating 
a conversion factor in XSPEC using the temperature in the 
[0.15-1] /?50() aperture, and subsequently corrected to take into 
account radial variations of temperature and abundance to give 
the final deprojected, PSF-corrected radial density profile. 

2.2.2. Temperature profiles 

Projected 2D temperature profiles were derived from spectra ex- 
tracted in logarithmically-spaced annular bins centred on the 
peak of the X-ray emission. Binning was such that the first bin 
was defined to have a significance of 30o" above background, and 
subsequent bins were defined so as to have RouilR'm - 1.33-1.5 
depending on the quality of the observation. The instrumental 
and particle background was subtracted from each annulus us- 
ing custom stacked, recast data files accumulated from obser- 
vations obtained with the filter wheel in the CLOSED position 
(FWC), renormalised using the count rate in a high energy band 
free of cluster emissior0. After subtraction of the FWC spectra, 
all spectra were grouped to a minimum of 25 counts per bin. 

We then fitted the FWC -subtracted spectrum of an annular 
region external to the cluster emission with a physically moti- 
vated model of the cosmic X-ray background, consisting of two 
MeKaL components pl us an absorbed power law with a fixed 
slope of r = 1.4 (see iLumb et alJ I200I lOe Luca & Molendil 
I2OO4I) . Spectra were fitted in the [0.3 - 10] keV range using 
statistics, excluding the [1.4-1.6] keV band (due to the Al line 
in all three detectors), and, in the EPN, the [7.45-9.0] keV band 
(due to the strong Cu line complex). In these fits the MeKaL 
models are unabsorbed and have solar abundances, and the tem- 
perature and normalisations are free parameters; the powerlaw 
component is absorbed by the Galactic absorption in the direc- 
tion of the cluster and since it has a fixed slope, only its normal- 
isation is an additional free parameter in the fit. This best fitting 
background model, with renormalisation appropriate to the ratio 
of the surface area of the extraction regions (corrected for chip 
gaps, etc.), is then added as an extra component in each annular 
fit. This is our standard fit. 

We then vary the power law index in such a way as to mimic 
a +10 per cent variation in the [2-10] keV flux and refit the spec- 
trum of the external region. The annular spectra were refitted 
using this new cosmic X-ray background model, and the change 
in the cluster temperature in each annulus relative to the standard 
fit was treated as a systematic uncertainty and added in quadra- 
ture to the statistical errors in each annulus. 

To deproject and PSF-correct, we assume that the 3D tem- 
perature profil e can be represented by a parametric model 
(adapted from IVikhlinin et al.l |2006|) that is convolved with a 
response matrix which simultaneously takes into account pro- 
jection and PSF redistribution. This model was projected, taking 
into acc ount the weighting sch eme proposed by Vikhlinin ( 200a, 
see also iMazzotta et alJl2()()4l) to correct for the bias introduced 
by fitting isothermal models to multi-temperature plasma emis- 
sion, and fitted to the observed 2D annular temperature profile. 
Uncertainties were estimated from Monte Carlo randomisation 
of the projected temperature profile assuming a Gaussian distri- 
bution defined by the uncertainties on each data point, and then 
corrected to take into account the fact that parametric models 
tend to over-constrain the 3D profile. Full details of the depro- 



^ Our adoption of FWC data allows the use of a physical model for 
the X-ray background, i n contrast to our previous analysis which used 
a blank sky background ( iPratt et aLll2007l) . The results are consistent. 
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Fig. 1. Entropy profiles of the R EXCESS sample, colour coded according to spectroscopic temperature measured in the [0.15-0.75] R500 aperture. 
Solid lines show the profiles derived from direct measurement; dotted lines show the entropy extrapolated into the central regions assuming an 
isothermal distribution at the temperature of the inner 3D data point (see text). On the left, the profiles are plotted in physical units; on the right, 
they are plotted in units of scaled radius R^qo estimated from the Msqq-Yx relation given in Eqn.[T] In the right hand plot, dashed lines mark, from 
left to right, radii corresponding to 0. 1 R200, ^250o(~ 0.45 R500) and Rioooi~ 0.7 Rjoo). 



jection and PSF correction of the temperature profiles, plus ex- 
tensive tests of the robustness of the method, will be detailed in 
a forthcoming paper 

2.2.3. Entropy profiles 

Since the density profiles are determined on a radial grid of sig- 
nificantly higher resolution than that of the temperature profiles, 
we determined the best fitting parametric 3D temperature profile 
on the same grid as that of the deprojected, deconvolved density 
profile and calculated the entropy, K - kT In^J^ , accordingly. 

In all cases, in the very central regions a single tempera- 
ture bin encompasses a region covered by several density pro- 
file bins (the median number is 5). Given that the central density 
of the galaxy cluster population exhibits a dispersion of up to 
two orders of magnitude and the overall density profile changes 
by up to th ree orders of m agnitude from the centre to the out- 
skirts (e.g.. I Croston et aLll2p08[). while the temperature varies 
only by a factor of 2-3 (e.g.. lPratt et al.ll2007h . it is clear that the 
characteristics of the density drive the properties of the entropy 
profiles. In order to examine the behaviour of the central en- 
tropy, we assume a constant central temperature, with the value 
given by the 3D temperature of the first shell. A simi lar pro- 
cedure was used by Donahue et al. (2005) and Cavagnol o et al.l 
(i2009t) and our adoption of this approach allows us to compare 
directly with their results. Note that for systems with poor central 
temperature profile resolution this extrapolation is only weakly 
model dependent, since it essentially concerns the disturbed sys- 
tems, whic h have rather flat c entral temperature profiles (see 
Figure 3 of lArnaud et allbOOgl) . When this scheme is applied, 
the XMM-Newton profiles have a typical central resolution of 



~ 5 /j^Q kpc, which compares favourably with the typical reso- 
lution of ~ 2 kpc in the Chandra analysis of Cavagnolo et al.l 
( l2009h . The left hand panel of Figure[T]shows these entropy pro- 
files plotted in physical units (/z^q kpc). 



2.3. Scaling 

In order to compare cluster profiles on a common radial scale, we 
express them in terms of Rg, the radius within which the mean 
mass density is 6 times the critical density at the cluster redshifQ. 
For practical purposes, we generally scale to ^500, the effective 
limiting radius for high quality observations from XMM-Newton 
and Chandra. Since the sample contains systems in a variety of 
dynamical states, we use Yx as a mass proxy. W e estimate /Jgoo 
iteratively as described in iKravtsov et alj (l2006i) . from the up- 
dated calibration of t he Ms qq-Yx relation obtained by combin- 
ing the lAmaud et al.l (l2007l) results from nearby relaxed clusters 
with REXCESS data from morp hologically relaxed s ystems. The 
full sample of 20 objects (8 from lArnaud et al.l2007l and 12 from 
REXCESS) is comprised of all systems for which the mass pro- 
files are measured at least down to a density contrast 6 - 550. 
The resulting M500 - Yx relation is: 



10 



14.567+0.010 



Yx 



2x 10"*MokeV 



iho Mo, 



(1) 



consistent with that derived by lArnaud et al.l (2001) b ut with 
impro ved accuracy on slope and normalisation ( Arnau det al.l 
i2009b . We also use the spectroscopic temperature T, measured in 
the [0. 15 - 0.75] R5Q0 region, to investigate the scaling properties 



Ms = 6pc{z)(47T/3)Rl, where p,(z) = E^{z)3Hl/8nG. 
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Fig. 2. Dimensionless entropy profiles of the R EXCESS sample compared to theoretical expectations from non-radiative simulations. The ob- 
served profiles have been renormalised by the characteristic A'joo as defined in Equation[3] Line styles are as for Figure[T] The dashed line depicts 
the best fitting power law fit to the the median entropy profile in the radial range [0.1 - l]i^200 for the clusters formed in the non-radiative sim- 
ulations of IVoit et al.l i2005l) . The dot-dashed line represents the same relation corrected for a 13 per cent underestimate of true mass due to the 
use of hydrostatic equilibrium. Left panel: clusters arranged according to temperature. The thick black line represents the median of all observed 
dimensionless profiles. The shaded grey area corresponds to the region enclosed by the median profile and typical scatter of the SPH simulations 
in lVoit et alj ( l2005h . Right panel: subsamples identified, defined as in Section [ZT] Black profiles denote clusters that are neither cool core nor 
morphologically disturbed. 



of the entropy and associated profiles. These values are given in 
Table [3 

The right hand panel of Figure [1] shows the entropy pro- 
files plotted in terms of R^qq. Plotting them this way explicitly 
shows the temperature dependence of the entropy distribution in 
the outer regions, and that despite the wide range of central en- 
tropy values, clusters clearly become more structurally similar 
with increasing radius. Beyond 0.2/?5oo, the relative dispersion 
in scaled entropy profiles is constant at approximately 30 per 
cent. 

In the pure gravitational collapse scenario, the scaled pro- 
files of any physical quantity should coincide, and so measures 
of these quantities at any scaled radius should correlate with 
global cluster parameters such as temperature or mass. The stan- 
dard self-similar scalings are A" oc Eiz)^'^^^ T for temperature and 
K oc £'(z)-2/3 for mass. Departures from these expectations 
are the direct result of the action of non-gravitational processes 
on the thermodynamics of the ICM. For comparison with previ- 
ous work on entropy scaling relations, the entropy was measured 
for all clusters at radii equivalent to 0.1 R200, ^2500 (~ 0.45 R5Q0) 
and Rim) (~ 0.72 Rsoi)), estimated from the scaling relations pre- 
sented in lArnaud et al.l (l2005h . via interpolation in the log-log 
plane. In addition, if R500 falls within the radial range encom- 
passed by the centre of the outer temperature annulus (13 sys- 
tems), we also calculated the entropy at Rsm. Uncertainties on 
the entropy were estimated from the quadratic sum of the errors 
associated with the deconvolved density and temperature pro- 
files. These values are listed in Table[T] 



3. Comparison to theoretical expectations 

3.1. Normalisation 

We first address the question of the absolute entropy normal- 
isation with respect to theoretical expectations in the case of 
purely gravitational structure formation. Any deviations from 
the predicted normalisation would point to the influence of 
non-gravitational processes on the thermodynamics of the ICM. 
Furthermore, radially resolved entropy profiles can be used to 
assess the radial extent of the change in ICM properties due to 
non-gravitational processes, if any. In the following, we compare 
the normalisation of our observed entropy profiles to theoretical 
expectations in order to elucidate the mass dependence and ra- 
dial extent of these effects. 

Numerical simulations which only implement gravitational 
processes make very specific predictions for both th e normal- 
isation and shape of galaxy cluster entropy profiles. IVoit et all 
(2005) discuss such simulations, pointing out that once scaled 
by the characteristic entropy of the halo. 



■'^200 - 



1 



27r G^M- 



200 



15 fbHiz) 



2/3 



(2) 



where fi, is the baryon fraction, the simulated SPH profiles, 
when fitted in the [0.1 - 1] R200 radial range, scatter about a me- 
dian scaled profile described by the baseline power law relation 
K(R)IK2(xi - 1.32 (/?//?20o)' ' with approximately 20 per cent 
dispersion. We note that non-radiative simulations do not pre- 
dict power law profiles down to arbitrarily small radii; in fact. 
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Fig. 3. Dimensionless entropy K/Ksoo as a function of mass M500 for different fractions of scaled radius. The solid lines show the power law fit to 
the data at each scaled radius with exponent a; the dashed lines give the expected values from the power law relation K/Ksqq = 1.42 (R/Rsoq)''^ 
fitted to the non-radiative numerical simulations of lVoit et alj ([2005). The band indicates the typical dispersion of the simulated profiles about this 
relation at the radii indicated. The dot-dashed lines show the same relation adjusted for a 13 per cent mass underestimate due to the assumption 
of hydrostatic equilibrium. Cool core and morphologically disturbed subsamples are indicated. There is a clear segregation in central entropy 
properties of the subsamples at 0.1 ^200. which disappears at larger scaled radius. The dimensionless entropy at R^oo is slightly higher than, but 
approximately consistent with expectations from simple gravitational collapse, independent of mass. 



significant flattening is generafly found interior to R < O.lRsoo- 
In addition, grid-based simulations consistently find larger en- 
tropy cores than SPH simulations even in the non-radiative case, 
which appears to be due to diff'erences in particle mixing be- 
tween the two diff'erent computational methods ( Mitchell e t aTl 
|2009|) . However, beyond ~ 0.1 R2 Q0, the agreement between sim- 
ulatio ns is better than 10 per cent (iVoit et al.l2005l:lMitcheU et all 
120091) . In the following, since we only compare at and beyond 
0.1 R20Q, we compare to the SPH results only. 

Assuming an abundance of Z = 0.3 Zq and mean molecular 
weight jj. = 0.596, the characteristic entropy can be rewritten: 



106 keV cm" 



10'* h-^ Me 



2/3 



2/3 

J,^ '''' 



-2/3 1,-4/3 



(3) 



For our assumed ft, — 0.15 and Rsm/Rioo - 0.659 derived from 
an NFW profile with csoo - 3.2 (the me an c measured for a 
morphologically relaxed cluster sample by iPointecouteau et all 
120051) . the baseline relation becomes 



K(R)/K5m^ 1 .42 (R/Rsoo) 



1.1 



(4) 



Note that this baseline relation was derived from simula- 
tions where the true masses were known, while our masses 
are calculated from an M-Yx relation derived from hy- 
drost atic mas s estima tes of mo rphologically relaxed sys- 
tems. Piffaretti & Valda rninil (l2008l and references therein) and 



lArnaud et alT(f2009). applying hydrostatic mass estimates to a 
large number of simulated clusters, argue that such masses can 
be underestimated by -13 + 16 per cent on average. In this case 
the normalisation factor in Eqn. |4]becomes 1.47. 
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Radius C a a]„K 

(keV cm"-) raw int 



Entropy-temperature relation 

0.1i?2oo 347 ±23 0.89 ±0.15 0.262 ± 0.040 0.254 ±0.041 

i?25oo 783 ± 15 0.76 ±0.06 0.120 ±0.021 0.083 ±0.118 

Jfiooo 1152 ±27 0.83 ±0.06 0.093 ±0.015 ... 

ifjoo 1489 ± 125 0.92 ± 0.24 0.265 ± 0.055 . . . 

Entropy-mass relation 

^2500 864 ±27 0.42 ±0.05 0.155 ±0.025 0.136 ±0.031 

i?iooo 1308 ±52 0.48 ±0.04 0.119 ±0.017 0.052 ±0.162 

«5no 1748 ±237 0.62 ±0.17 0.265 ± 0.055 0.221 ±0.160 



T is the spectroscopic temperature in the [0.15 - 0.75]i?joo region; masses are estimated from the M - Yx relation given in Eqn.[T] Data were 
fitted with a power law of the form E(z)"K = C x (A/Aq)", with Aq = 5 keV and 5.3 x 10'"* Mq, and n = 4/3 and 2/3, for temperature and mass 
repectively. Fits used orthogonal BCES regression with errors estimated using bootstrap resampling. The raw and intrinsic logarithmic scatter 
about the best fitting relations are given in the final two columns. 



Figure |2] shows the dimensionless entropy profiles 
K{R)/K5Qo- Clearly, the observed profiles do not coincide. 
The central regions show the most dispersion and the profiles 
tend to converge towards the non-radiative prediction at large 
radius, but those of the lowest mass systems converge slowest; 
in other words, their slopes are shallower (discussed in more 
detail below). This is a manifestation of the fact that the entropy 
modification extends to larger radii in lower temperature 
systems, consistent with the expectation that non-gravitational 
processes have a greater effect at the low mass end of the cluster 
population. 

To better quantify the above, Figure[3]shows the dimension- 
less entropy K/Ksqi) versus mass for various fractions of scaled 
radius. Also overplotted in each panel is the expectation from 
Equation |4] As expected, it can be seen that at 0.1 Rsqo, the ex- 
cess with respect to the theoretical prediction from gravitational 
collapse is strongly mass-dependent, with the least massive sys- 
tems exhibiting the strongest deviation. The mass dependence 
becomes less pronounced as we proceed towards the outer re- 
gions of the ICM, such that at R5Q0, the mass dependence is en- 
tirely consistent with zero and the measured values scatter about 
the theoretical prediction. Indeed, at this radius, the very slight 
negative slope can be attributed to the single lowest mass data 
point, which drives the fit. 

At small radii there is large amount of scatter ((T\nKM - 
0.26 + 0.04 about the best fitting regression line). Dividing the 
data into subsamples elucidates the origin of this scatter; there is 
a clear segregation in subsamples, with cool core systems show- 
ing the least deviation from the baseline prediction while mor- 
phologically disturbed systems show the most deviation. The 
subsample segregation disappears as we push outward though, 
and there is no evidence for any segregation at or beyond Rum- 
The full radial behaviour of this trend is explicitly illustrated in 
the right hand panel of Figure |2l 

Interior to /?5oo, the observed entropy is always higher than 
the baseline prediction. However, at Rsa), the median dimension- 
less entropy is /r(7?5oo)/^500 = 1.70+0.35, where the uncertainty 
comes from the standard deviation of the points. This is slightly 
higher than, but consistent with both the baseline prediction of 
KiRsoi)) / K^iK) = 1.42 (Eqn.lU and the same prediction corrected 
for a 13 per cent mass bias due to the assumption of hydrostatic 



equilibrium, K(Rs()q) / Ksqq = 1.47. The lack of mass dependence 
and agreement with the normalisation from simulations was als o 
noted for a sample of cool core clusters bv iNagai et al.l ( l2007h . 
For the present representative sample, given the large uncertain- 
ties it is not possible to test the predictions more thoroughly, 
underlining the need for robust, high quality, spatially resolved 
entropy measurements at and beyond /?50(). 



3.2. Entropy scaling relations 

For comparison with previous work, we also examined the 
entropy-temperature and entropy-mass relations. Scaling rela- 
tions were fitted with a power law of the form E{z)"B = 
CiA/Ao)", with Ao = 5 keV and 5.3 x 10^'^ for T and M 
respectively, and n fixed to the expected scaling with redshift 
(« - 4/3 for T and 2/3 for M). Data were fitted u sing the orthog- 
onal BCES minimisation technique dAkritas &^B ershadv 199g) 
with uncertainties on each fit parameter estimated from bootstrap 
resampling. 

The best fitting slopes and intercepts for the entropy- 
temperature and entropy-mass relations at various scaled radii 
are listed in Table |2l The evolution of these slopes with increas- 
ing radius mirrors the behaviour of the dimensionless entropy 
discussed above; in the inner regions the relations are shallower 
than self-similar with large scatter, while at R^cm the relations 
are compatible with self-similar (although with relatively large 
uncertainties given the limited number of data points). 

Comparing to previous work, a wide variety of slopes have 
been found from fits to the entropy-temperature relation at 
0.1 ^200^ ranging from ve ry shallow (a - 0. 49 + 0.15; |Pratt et al.l 
20061; a = 0.50 + 0.08: [Nagai et alj |2007l) to very steep (a = 
0.92 + 0. 12: [Sanderson et aI.ll200Sr g = 0.85 + 0.19: this work). 
We simply note that cool core-only samples tend to yield shal- 
lower slopes than statistically-selected samples, a fact borne 
out in the present data, for which the entropy temperature re- 
lation at 0. 1 R200 for the cool core subsample has a slope of 
a - 0.63 ± 0.94, while the morphologically disturbed subsample 
has a slope of ff = 1.22 + 0.76. Beyond 0.1 /?2()0, both the slope 
and the normalisation of the relations are in very good agreemen t 
with recent determinations ( Nagai et al.ll2007l; [Sun et alj|2009l) . 
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showing the excellent consistency between XMM-Newton and 
Chandra result^ 

None of the studies listed above give constraints on entropy 
evolution. In this context, our results underline the need for 
representative samples, to establish the effects of non gravita- 
tional processes and dynamical state on the evolution in the cen- 
tral regions. Furthermore, precise measurements at large radius 
(R > /?250o) are needed to establish the baseline entropy evolu- 
tion in the absence of non-gravitational effects. 



4. Radial entropy structure 

Various semi-analytical models and cosmological simulations of 
clusters formed in the absence of non-gravitational processes 
have shown that outside the central regions {R > 0.1 R200, or ~ 
0.15 j^50o) , entropy profi l es follow a power law with K{R) qc j^' ' 
dTozzi & NormanI I2 OO II: iBorgani et al.1 l2005t IVoit et alj llool 



Table 3. Entropy profile parameterisation with a model of the form 
K{R) = Kq + A'loo (^/100/!:^(J kpc)". Entropy is measured in units of 
keV cm"^. 



I2OO5I: iMitchell et a l. 2009). Simulated p rofiles flatten in the 



very cen tral regions due to entropy mixing dWadsley et al.ll2068t 
[Mitchell et al.ll200"9h . Observed profi les are also found generally 
to ha ve similar external slopes (e.g.. iPratt et al.ll2b06l; ISun et al.l 
l2009i) and to flatten in the central regions in high resolution 
Chandra observations dPonahue et all 120061: [Sanderson et al.l 
I2009t ICavagnolo et al.ll2009 r In the following, we investigate 
the central entropy slope and a parameterisation of the entropy 
profiles, and relate this to global cluster properties. 



4.1. Central slope 

The left hand panel of Figure |4] shows the probability density 
distribution of the logarithmic slope dlnK{R)/dlnR of the en- 
tropy profiles measured at 0.075 /?5oo, where we have directly 
measured data for all but two systems. The distribution shows 
two peaks, with the larger peak at a slope of ~ 0.4, containing 
approximately 2/3 of the sample, and a smaller peak at a slope 
of ~ 0.9, which comprises the cool core systems defined in the 
subsample classification scheme discussed in Section 12.11 Two 
peaks in the distribution of inn er logarithmic slope w ere also 
found in the Chandra analysis of I Sanderson et alj (l2009l in their 
case, at 0.05 Rsqo)- 



4.2. Parameterised fitting 
4.2.1. Model 

A power law is the simplest parameterisation of the entropy pro- 
files, but Figure [T] shows that such a description is clearly in- 
adequate to describe the majority of the profiles in the present 
sample if the core regions are included in the fit. A possibly 
more interesting parameterisation is to instead model the en- 
tropy profiles with a power law plus a c onstant, K(R) = Ko + 
Km{R/lQQh-^kpc)" (as introduced by iDonahue et all l2005l) . 
where Kq represents the typical excess of core entropy above the 
best fitting power law at large radii. In this case, to ensure that 
the same radial range is fitted in all cases, and to ensure the pre- 
cision of the powerlaw fit at large radii, the profiles were fitted 
between - ^1000] ■ The results of fitting the power law plus 
constant model are listed in Table [3] The external power law is 
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^ Our earlier results suggested somewhat shallower relations at i^iooo 
dPratt et alj|2006i . The difference can be traced to the increased preci- 
sion on the temperature profiles afforded by the present data, especially 
at low mass. 



extremely stable to fitting ranges and clearly does not depend on 
the entropy of the central region or the value of A^ioo- 



4.2.2. Distribution of central entropy, Kq 

The histogram of Kq values for the REXCESS sample is shown 
as a probability density in the right hand panel of Figure|4] where 
the bins are 0.15 dex and clusters with a Kq consistent with zero 
(i.e., pure power law profiles) are shown at the extreme left of 
the plot. For the REXCESS sample, the number of systems with 
a Kq consistent with zero at 3cr significance is three for a con- 
stant centalJemger^ireJlO per cent of the sample), consistent 
with ICavagnolo et al.l (l2009l) . who analysed a large number of 
Chandra archive observations and who also used a constant cen- 
tral temperature assumption. 

ICavagnolo et al.i (2009) found that the distribution of central 
entropies in their sample is bimodal with peaks of approximately 
the same amplitude at Kq ~ 15 keV cm^ and Kq ~ 150 keV 
cm^, and a distinct gap between /To ~ 30 - 50 keV cm^ (see 
also Rei prich & Hudsonll2007b . The peak values are indicated in 
Figure |4] by a dashed line. The distribution of REXCESS cen- 
tral entropies also exhibits two peaks, although there are subtle 
differences in their positions and amplitudes. In particular, there 
are more clusters with a high Ko than with a low Kq, and the 
positions of both peaks are shifted somew hat to lower Kq with 
respect to those found bv ICavagnolo et all 
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Fig. 4. Central regions of the REXCESS entropy profiles. Left panel: Probability density plot of the logarithmic slope of the density profile 
measured at a fiducial radius of 0.075 Rsoo- Clusters with a steep slope correspond to cool core systems. The solid line is a kernel density plot with 
a smoothing width of 0.15. Right panel: Probability density plot of the central entropy excess above a power law, Kq, from a power law plus a 
constant model fit K(R) = Kq + Kiqq (R/lOOh^g kpc)", to the entropy profiles. The solid line is a kernel density plot with a smoothing width of 0.1. 
Cool core clusters have the smallest values of Kn. 
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Fig. 5. Outer regions of the REXCESS entropy profiles. Left panel: Probability density plot of outer slope values, a, from a power law plus a 
constant model fit K(R) = Kq + Kioo(R/lOOhjQ kpc)", to the entropy profiles. The dashed line shows a value of 1.1. The solid line is a kernel 
density plot with a smoothing width of 0.1. Right panel: Entropy profile slope vs temperature in the [0.15 - 0.75]i?5oo aperture. Cool core and 
morphologically disturbed subsamples are indicated. 



4.2.3. Entropy slope outside the core 

The left hand panel of Figure |5] shows the histogram of fitted 
values of the outer profile slope a obtained from the power law 
plus constant model. While there is a quite substantial spread of 
values in entropy profile slope, ranging from extremely shallow 
(a ~ 0.5) to extremely steep (a ~ 1 .9), there is no indication for 
bimodality in the distribution of outer entropy slope when the 
profiles are modelled in this way. The median slope for the power 
law plus constant model fits is 0.98, which is slightly shallower 
than the canonical value of 1 . 1 (indicated by the dashed line in 
th e Figure). I t is low er still than the value of 1 .2 which was found 
by I Volt et all (|2005|) when fitting the [0.2 - 1] radial range. 



The right hand panel of Figure |5] shows the distribution of 
slopes versus temperature. The Spearman rank correlation co- 
efficient is 0.53 with a significance of 2.2 x lO"-' indicating a 
significant correlation of slope with temperature. This is most 
likely a manifestation of the well-known dependence of outer 
density profile slope on temperature (see Croston et al. 2008, for 
the present sample). The different subsamples are indicated in 
the Figure by blue stars (cool core systems) and red squares 
(morphologically disturbed systems). Clearly the cool core sub- 
sample has a very small scatter in outer entropy slope; in addition 
there is a more pronounced trend with temperature, illustrated by 
the fact that the Spearman rank coefficient for this subsample is 
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0.77 with a significance of 9.7 x 10 The morphologically dis- 
turbed subsample has a very large scatter, incorporating both the 
upper and lower extremes in outer slope values. The tight range 
of slopes in the cool core systems, together with the wide range 
of slopes in the morphologically disturbed systems, suggests that 
the core and external properties of the ICM are linked. 

5. Discussion 

5.1. Bimodality? 

Two peaks are seen in the central regions of the R EXCESS sam- 
ple, which are visible both in the distribution of logarithmic en- 
tropy slope at 0.075/?50() and in the distribution of Ko from pa- 
rameterisation of the profiles with a power law plus constant 
model (Fi gure m. For a more qu antitative comparison with the 
results of ICavagnolo et all ( 12009 ). we performed a maximum 
likelihood fit of the 29 clusters for which Kq is constrained. The 
fitting was performed on the unbinned data in log space, using 
the MCLUST and fitdistr packages in version 2.9 of the r statisti- 
cal software environment. We considered three different mod- 
els: a single Gaussian, a left-skewed distribution, and a double 
Gaussian (i.e., a bimodal distribution), and used the Bayesian in- 
formation criterion (BIC. ISchwarzll 1 9781) to distinguish between 
them. A difference in BIC of between 2 and 6 indicates posi- 
tive evidence against the model with the greater BIC value; val- 
ues above 6 indicate strong evidence against the model with the 
greater BIC value. The BIC values ai-e -67.53, -63.32 and -63.08 
for the single, left-skewed and double Gaussian distributions, re- 
spectively. Thus while the single Gaussian distribution is clearly 
the worst description of the data, on the basis of this test, and 
given the limited number of data points at our disposal, we can- 
not definitively distinguish between a bimodal and a left-skewed 
distribution of Kq. 

The best bimodal distribution fit yields Gaussian means at 
Kq ~ 3 and Kq ~ 75 keV cm^, with an amplitude ratio 
of 1:3. If we as sociate these values with the peaks found by 
ICavagnolo et al.l ( |2009) , then they are somewh at offset to lower 
values and the amplitude ratio is different (iCavagnolo et alj 
found peaks at ~ 15 and ~ 150 keV cm^, with an amplitude ratio 
close to 1:1). The shift of the higher central entropy peak is due 
to a number of clusters in the R EXCESS sampl e that fall direct ly 
in the Kq ~ 30 - 50 keV cm^ gap found by Cavagnolo et al.l It 
may be that these clusters have a more 'typical' value of Kg for 
the general population. Given that the ICavagnolo et al.l (l2009l) 
sample was archive limited, and given that prevailing sociologi- 
cal trends in cluster research have for many years led to a focus 
on extreme cool cores and spectacular mergers to the exclusion 
of more mundane objects, it may well be that clusters with a 
more 'typical' central entropy do not exist there but do exist in 
a representative sample. We note however that Cavagnol o et alj 
also detected bimodality in a complete flux-limited subsample 
of their data (although at much reduced significance), with the 
peaks shifted slightly lower relative to those of the full archive 
sample. Alternatively, the recent 2009 January Chandra calibra- 
tion updat e may offer a mo re prosaic, albeit partial, explanation. 
Since the ICavagnolo et al.l analysis predates this update, hotter 
systems with flat central temperature distributions {kT > 4 keV), 
would have a systematically higher temperature, and hence en- 
tropy. 

The shift of the lower central entropy peak appears to be due 
to technical differences con nected to the treatment of tempera- 
ture profiles. 'Cav agnolo et al.i derived their entropy profiles from 
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Chandra data, which has sub-arcsecond resolution in the cen- 
tre of the detector, meaning that no PSF correction was needed. 
However, while their density profiles were deprojected, their 
temperature profiles were not. For flat central temperature distri- 
butions, this will not substantially change the resulting entropy. 
However, for steeply declining central temperature distributions, 
neglect of projection effects will tend to lead to an overesti- 
m ate of the temperatu re of the inner annulus (see e.g.. Fig. 8 
of lPointecouteau et al.l l2004). and thus of the central entropy for 
the assumption of constant central temperature. The net effect 
of neglecting deprojection would be to shift the peak in Kq to 
higher values in the Chandra analysis, as observed. 

We note that the d ifferences in the central entropy distribu- 
tions of R EXCESS and lCavagnolo et al.l(l2009l) are unlikely to be 
due to resolution effects. We extracted surface brightness profiles 
in 3'.'3 bins while Cavagnolo et al. (2009) used 5" bins, meaning 
that for a given redshift, our XMM-Newton entropy profiles in 
fact have a higher resolution than those from Chandra. Our 3'.'3 
bins give us a resolution of 3.5-10 h^^ kpc for the present sam- 
ple, depending on redshift, with a median value of 7 /i^q kpc. 
Given that entropy cores extend ty pically to 30 kpc in C handra 
observations of cool core systems (iDonahue et al.ll2005h . these 
should be easily detectable with the current data if they exist. 
However, Chandra follow-up observations would still be desir- 
able to quantify the cores of these systems at higher resolution, 
and particularly to investigate if there is indeed a turnover in en- 
tropy at very small radius in systems where we are unable to 
detect this effect. 

We also do not find any bimodality in the distribution of 
outer power law slope a, consistent with the increased self- 
similarity of the profiles in the outer regions. In this context it 
is important to underline the fact that our fitted model consisted 
of a power law plus an additive constant, and is thus different 
from t he simple power law model fitted by e.g.. I Sanderson et aTl 
(l2009l) . 



5.2. Maintaining the distribution of centrai entropies 

The clear link between the lack of a cool core and the presence of 
morphological disturbance established by the R EXCESS sample 
gives important clues to the processes at play in giving rise to the 
observed distribution of central entropies. It seems paradoxical 
that the profiles of cool core systems, where non-gravitational 
processes play an important role, resemble most closely the non- 
radiative baseline, while those of the unrelaxed objects, whose 
properties are expected to be dominated by gravitational pro- 
cesses, deviate the most from the baseline. 

A possible explanation for the form of the cool core profiles 
is that in these systems, AGN heating is gentle and serves pri- 
marily to balance cooling, thus preserving the increasing form 
of the entropy profiles. If this were indeed the case, it would 
imply that AGN heating is achieved via mechanisms such as 
weak shocks (e.g., Fabian et al . 2003), or buoyant bubbles (e.g., 
iChurazov et al.ll200ll) . rather than via catastrophic explosions. 

We can envisage two scenarios that could explain the proper- 
ties of the morphologically disturbed clusters in R EXCESS and 
the link between the lack of a cool core and the evidence that 
a cluster is unrelaxed. One explanation is that a combination of 
extra heating and dynamical activity have conspired to keep the 
central entropy elevated and so prevent the morphologically dis- 
turbed systems in the R EXCESS from ever forming a cool core 
in the first place. This is possible because post merger distur- 
bance persists for longer in these clusters, as in a disrupted core 
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Fig. 6. Implied Spitzer conductivity suppression factor as a function of 
radius. The different cluster subsamples are indicated; line styles are as 
for Figure[T] The dashed line indicates the threshold between thermally 
unstable and conductively stable regimes. 



Fig. 7. Gas mass fraction profiles /gas(< R)- Total mass profiles 
are calculated for an NFW profile with a concentration C500 = 3.2 
dPointecouteau et al] l2005f) , normalised to the M500 estimated from 



with less entropy contrast, buoyancy differences will not be as 
strong, and the stratification that restores a cluster to a relaxed 
state will happen more slowly. 

The second explanation is simply that the morphologically 
disturbed systems in the R EXCESS were originally cool core 
systems, and that mergers have disrupted the entropy structure 
in the central regions. Given time, they should relax back to 
a cool core state, but we are seeing them before the cool core 
can re-establish itself. We note that such a scenario is supported 
in terms of timescales by the observed relationship between the 
lack of cool cores, clear evidence of dyna mical disturbance an d 
the presence of radio halo emission (e.g. iGovoni et alj|2004l) . 
The ongoing REXCESS radio follow-up will be extremely use- 
ful in determining if this possibility can explain the clear anti- 
correlation of cool cores and disturbed morphology. 

Furthermore, it is possible that the timescale for cooling 
to (re-) establish itself is impaired by conduction. Figure |6] 
shows the conduction suppression factor fc expressed in terms 
of the radius (in kpc) and the entropy (in keV cm^), viz., fc ~ 
62.5 {R^IK^) (Sanderson et al. 2009ll|. It is clear that non-cool- 
core clusters can be stabilised since conductivity would need to 
be suppressed at most by only a factor of fc ~ 0.1 to locally 
counteract cooling (Figure while this will not be effective 
in the centres of most cool core systems, where the threshold 
for conductive stability (fc - 1) is e xceeded (a similar result 
was found by ISanderson et al.l l2009l) . Thus conduction could 
also contribute to sustaining the elevated central entropy in post- 
merger systems long after the initial disruption. 



' Note that this equation only holds at the limit of conductive thermal 
balance, i.e. where the radius exactly matches the Field length at all 
radii. 



We note that high central entropy non-cool-core systems are 
difficult to reproduce in numerical simulations, which seems to 
be a consequence of the presence of very dense cores which are 
very hard to disrupt (e.g. Gomez et al. 2002; Poole et al. 2008j). 
The presence of these dense cores, of which in some cases mul- 
tiple instances may be present in the same system (Bohringer et 
al., in preparation), at variance with observations, may point to 
deficiencies in the modelling of the complex interplay between 
gravitational and non-gravitational processes in these simula- 
tions. 

5.3. Linking entropy and gas mass fraction 

Entropy modification is generally discussed in terms of three 
basic mechanisms: early heating ('pre-heating'), where the 
gas is heated before accretion into the dark matter potential 
well, pres umably either by early supernovae a nd/or AGN ac- 
tivity (e.g. lKaise3ll991l: lEvrard & Henirvlll991h : internal heat- 
ing after accr etion by the same or similar mechanisms (e.g. 
Met zler & Evr ard 1994; Bower et al. 2008); radiative cooling of 
the gas (e.g. Pearce et al.l 120 00). where the lowest entropy gas 
found in the centre of the cluster condenses and cools out of the 
ICM. All of these processes act to change the total amount of gas 
in the central regions of a cluster, either through making it more 
difficult to compress into the halo (early heating), through con- 
vection of gas to the outer regions (internal heating), or through 
physical removal of the gas to form stars (cooling). 

In Figure I2] we show the gas mass fraction profiles /gas (< 
R) - Mgas (< R)/M{< R) for the present sampl e. Gas masses 
have been calculated from the gas density profiles dCroston et al.l 
2008). Total mass profiles were calculated assuming an NFW 
profile with concentration C500 = 3.2, the average concentration 
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Fig. 8. The dimensionless entropy K/Ksqq rises with fg(R) = 
Mg/ifbMsoo), the fraction of a cluster's baryons in the ICM within ra- 
dius R in all clusters. However, the entropy distribution shows a clear 
temperature dependence. The dashed line illustrates the unmodified en- 
tropy distribution for a cluster of M500 = 8 x 10'* Mq (approximately 
the mass of the most massive cluster in the present sample), assuming a 
concentration C500 = 3.2 and gas in hydrostatic equilibrium in the clus- 
ter potential with a density profile identical to that of the dark matter. 



derived from t he total mass profiles of the morphologically reg- 
ular sample of 'Po intecouteau et al.l (|2005|FI normalised to M500 
estimated from Eqn. [l] There is a clear dependence of /gas on 
temperature/mass, throughout the observed temperature range, 
in the sense that hotter, more massive systems have higher gas 
mass fractions throughout the ICM. In addition, there is a clear 
dependence of gas mass fraction with radius in all systems, and 
only the most massive clusters have gas mass fractions which 
approach the universal value at the highest radii we are able to 
probe. 

In Figure [8] we plot the dimensionless entropy K/ K^oq as a 
function of fg (< R) = Mg^s/ifh M^qq), the fraction of a cluster's 
baryons in the ICM within radius R. fg is calculated assuming 
fb = 0.15 (flfe/z^ - 0.022 and n„, - 0.3), and using total masses 
estimated from Eqn.[T] Overplotted for comparison is an unmod- 
ified entropy distribution from the models of Voit et al. (2002), 
derived for a cluster of M500 = 8 x lO''* Mq (approximately the 
mass of the most massive cluster in the present sample), assum- 
ing a concentration C500 = 3 and gas in hydrostatic equilibrium in 
the cluster potential with an identical density profile to that of the 
dark matter. This particular representation makes explicit both 
the dependence of the entropy distribution on baryon (gas) frac- 
tion, and the mass/temperature dependence of the baryon (gas) 
fraction itself. The gradual translation of the profiles to the left 



* The dependence of concent ration on total mass is negligible for 
the m ass range we consider here jPointecouteau et alj2005hlBuote et al] 
l2007b . 



Fig. 9. Entropy profiles multiplied by the gas mass fraction profile. The 
dashed line is the predicted entropy distribution from the non-radiative 
simulations of IVoit et al] ilOOS) . 



hand side of the plot is due to a systematic lack of baryons (gas) 
in low temperature systems relative to high temperature systems. 

The consequence of entropy modification is thus to remove 
gas (or prevent gas from accreting) in such a way as to leave 
both a radial and a mass dependence in the gas mass fraction. It 
is interesting to see whether connecting for this effect might bring 
the entropy profiles of our sample into agreement with the pre- 
dictions from simulations. In Figure |9] we show the dimension- 
less entropy profiles multiplied by the gas mass fraction profile 
/gas (< R), a renormalisation that is equivalent to correcting si- 
multaneously for the global and radial dependence of the gas 
mass fraction. Once renormalised in this way, the profiles are al- 
most fully consistent, both in slope and normalisation, with the 
expectation from non-radiative simulations, and the dispersion 
drops dramatically. Slightly better agreement can be found if the 
simulated relation is multiplied by a factor to take into account 
the hydrostatic mass offset found in numerical simulations. 

The above results can be explained in a number of ways. 
There may be a bias of gas accretion with mass, for instance due 
to early heating, which makes it more difficult to compress the 
gas into lower mass haloes. Once accreted, the gas may be re- 
moved from the hot phase by radiative cooling, which acts more 
efficiently in the densest central regions. However for this mech- 
anism to be wholly responsible for the observed entropy proper- 
ties would require it to affect the gas to a considerable fraction 
of /?5oo and to be preferentially efficient in low mass systems. 
Finally, there may be expulsion of material from the central re- 
gions towards the outskirts, perhaps via convection due to en- 
ergy injection from supernovae or a central AGN, with the low- 
est mass systems experiencing the most central gas removal due 
to their shallower potential wells. 

The mass dependence of the total baryon fraction provides a 
way of discriminating between the competing processes. If cool- 
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ing is the dominant effect, then the total baryon fraction should 
be almost constant across the mass range from groups to clus- 
ters as the low entropy gas is converted into stars. However, sim- 
ply cooling out the low entropy gas would greatly exceed the 
observed mass in stars, and would lead to a galaxy luminosity 
function completely at odds with observations. If instead extra 
heating dominates, then the baryon fraction should be lower for 
lower mass, group-scale systems, as either early heating makes 
compression of gas into these haloes less efficient or AGN activ- 
ity expels gas from their shallow potential wells. The observed 
anti-correlation in the relative dependencies of gas and stars with 
total mass implies that there is more ma ss locke d in stars in 
systems which contain les s gas (ILin et al.l i2003; Gonza lez et al.l 
120071: iGiodini et al.ll2009l) . However, recent results suggest that 
there is still a ~ 3cr deficit of baryons wit h respect to that mea - 
sured by WMAP on galaxy group scales (iGiodini et al.ll2009l) . 
implying that both cooling and heating must contribute to chang- 
ing the thermodynamic properties of the ICM. 



5.4. Speculative scenario 

The representative nature of the R EXCESS sample has brought 
to light some intriguing points outlined above, which allow us 
to propose a tentative scenario. It seems that about two thirds 
of the R EXCESS clusters possess a significantly higher central 
entropy than that expected from cuiTent non-radiative cosmolog- 
ical simulations and consequently do not possess a cooling core. 
A combination of extra heating and continuous ICM mixing due 
to merging may have kept these systems on a higher adiabat, 
leading to the observed high central entropies. Some early ex- 
tra heating may have occurred in the protocluster phase, which 
would coincide with the peak in AGN activity at z ~ 2 - 3. In 
this scenario the lower entropy envelope traced by the non-cool 
core systems (see Figure|2]) could indicate the level of early extra 
heating. The distribution of central entropy above this lower en- 
velope would then be produced by later heating and gas mixing 
during mergers, with the least relaxed objects having the higher 
central entropy, as observed. These processes will inhibit forma- 
tion of a stable cool core and naturally lead to redistribution of 
the gas to the outskirts, acting most importantly in low mass sys- 
tems, corresponding to the observed behaviour of the gas mass 
fraction. 

In contrast, about one third of the R EXCESS sample possess 
a cool core. The clear association of the BCG with the bottom of 
the potential well ( IHaarsma et al.ll200 9) and their regular X-ray 
morphology testifies to the relaxed nature of these objects. These 
systems presumably experienced a less chaotic early dynamical 
history leading to a modest entropy elevation due to mixing (if 
any), and may have undergone less early extra heating, allowing 
them to develop a cool core at a relatively young age. The natural 
reduction of entropy due to cooling while the gas is still in the 
hot phase, due to the combination of a temperature drop and 
the consequent increase in gas density needed to keep pressure 
balance, may bring the profiles into line with the observed power 
law behavior. However, cooling must be limited by some finely- 
tuned feedback process to prevent a significant fraction of the 
gas from disappearing from the hot phase, which would lead to 
a net increase in entropy, at variance with the observations. 

Our scenario bears some resemblance to the model proposed 
by McCarthy et al. (2008), although with increased emphasis on 
merger mixing as a process for setting and maintaining entropy 
levels in non-cool core systems. 



6. Conclusions 

Our data represent a considerable advance over those used in 
most previous analyses of the entropy structure and scaling in 
clusters. The sample of 31 clusters spans the temperature range 
[2-9] keV and includes systems with a variety of entropy char- 
acteristics. The objects have all been observed to approximately 
the same depth with the same instrument, allowing us to probe 
the properties of the entropy out to significant fractions of R^oo 
(at least /Jiooo for all systems, and at least out to Rsm for thirteen 
systems), which is essential to determine the radial extent of the 
effect of non-gravitational processes on the ICM. 

In the inner regions, there is a mass dependent entropy excess 
with respect to theoretical expectations derived from cosmologi- 
cal numerical simulations including only gravity. At larger radii, 
the mass dependence weakens and the dispersion drops dra- 
matically. The mass dependence disappears at R500, and the en- 
tropy normalisation is, within the relatively large observational 
and theoretical uncertainties, in agreement with the expectations 
from non-radiative numerical simulations. This behaviour is mir- 
rored in the entropy scaling relations, which are non-self similar 
at small radii but are compatible with self-similar at TJgoo. While 
si milar results were f ound for a sample of cool core clusters 
by iNagai et al] (|2007|) . it is important to note that in the group 
regime, at temperatures lower than 2 keV, the entropy normali- 
sation of morphologically relaxed s ystems has been found to be 
significantly higher than predicted (fSun et alj|2009l) . 

In the inner regions there is considerable dispersion, with 
a distinct segregation in residuals, in the sense that cool core 
clusters show the least deviation and morphologically disturbed 
systems show the most deviation with respect to expectations 
from non-radiative simulations. This dependence disappears at 
~ R\oQo- This clear association of unrelaxed morphology and el- 
evated central entropy would suggest either that cool cores are 
destroyed by mergers, or that cool cores have never been able to 
form in these systems. 

Fitting the entropy profiles with a power law plus constant 
model allows us to constrain a, the power law slope at large 
radius, and also Kq, the central excess of entropy with respect 
to this power law. With the current data we cannot statistically 
distinguish between a bimodal distribution or a left-skewed dis- 
tribution of Kq in log space; however, there is certainly no ev- 
idence for strong bimodality in the present sample. The distri- 
bution of outer entropy slopes is unimodal, with a median slope 
of 0.98. Cool cores have a narrow range of outer entropy slopes 
(0.8-1.2) while morphologically disturbed systems have a much 
wider range of outer slopes (0.5-1.9), suggesting a link between 
the properties of the cores and the outer regions of clusters. 

In seeking to explain the structural and scaling behaviour 
of the entropy we looked at the gas mass fraction profiles of 
the sample. These are strongly mass dependent. Furthermore, 
the gas mass fraction increases with radius in all cases. A plot 
of the dimensionless entropy versus baryon fraction explicitly 
shows the dependence of entropy on gas mass. Renormalising 
the dimensionless entropy profiles by the gas mass fraction pro- 
file /gas(< R), effectively correcting simultaneously for both the 
mass and radial dependence of the gas mass fraction, dramati- 
cally decreases the dispersion in scaled profiles and brings them 
into agreement with predictions from non-radiative simulations. 
This provides further evidence for the underlying regularity of 
the cluster population, which has important implications for their 
use as cosmological probes. 

The implication is that variations of gas content with mass 
and radius can explain the observed properties of the entropy 
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distributions, and by implic ation, the suppres sion of luminosity 
in low mass systems (e.g.. iPratt et al.ll2009h . Various physical 
mechanisms can impose this behaviour, but the mass and radial 
dependence would strongly argue for a combination of both ex- 
tra energy input and radiative cooling. However, it is necessary 
for these mechanisms to be capable of affecting the physical 
properties of the gas at least up to Rum, and perhaps beyond, 
in order to explain the observed offsets with respect to expec- 
tations from non-radiative numerical simulations. We discuss a 
tentative scenario to explain the observed behaviour of the en- 
tropy and gas mass fraction in the R EXCESS sample, based on 
a combination of extra heating and merger mixing maintaining 
elevated central entropy levels in the majority of the population, 
with a smaller fraction of systems able to develop a cool core. 

In the near future, our ongoing R EXCESS radio follow-up 
will help shed light on the relationship between merging activity 
and the lack of cool cores. Looking further ahead, observations 
of a similarly selected sample of galaxy groups, which are more 
sensitive to the effects of non-gravitational processes, would help 
to establish the magnitude of their impact at low masses. Precise, 
spatially-resolved measurements of the entropy at large radius 
(R > Rsoo) will be a further essential test of theoretical mod- 
els. Furthermore, a full census of the baryonic matter across the 
entire mass range for a representative sample of groups and clus- 
ters is essential to determine the relative contribution of the dif- 
ferent processes in play. 
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